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Recent  advances  in  communications  technology  suggest  that  more  powerful  systematic  methods 
for  the  reduction  of  IMG  must  be  developed  to  keep  pace  with  this  technology.  In  this  chapter  analyti¬ 
cal  approaches  to  IMG  diagnosis  are  reviewed  and  discussed.  It  is  found  that  a  great  deal  of  informa¬ 
tion  about  the  nature  of  the  IMG-producing  nonlinearity  can  be  deduced  from  the  intermodulation 
(intermod)  spectrum  and  this  can  provide  important  guidelines  in  identifying  the  source  of  the  non¬ 
linearity.  On  the  other  hand  it  is  not  possible  to  identify  with  certainty  a  unique  physical  source  of  an 
intermod  from  the  intermod  spectrum  alone  although  the  reverse  is  possible.  The  strengths  and  limita¬ 
tions  of  the  analytical  approach  are  delineated  and  the  need  for  the  application  of  high  speed  computers 
for  Fourier  transform  analysis  of  signals  containing  intermods  is  demonstrated.  Several  model  com¬ 
puter  calculations  are  carried  out  and  the  results  are  discussed  and  compared  when  appropriate  to  what 
might  be  expected  from  analytical  considerations.  Insight  is  provided  into  such  effects  as  the  power  of 
the  intermod  relative  to  primary  signal  power,  the  surprising  decrease  in  power  of  certain  intermods 
when  additional  primary  signals  are  turned  on,  the  relative  power  of  different  intermods,  etc.  The  goal 
is  to  provide  a  mathematical  and  computational  basis  for  predicting  the  intermodulation  spectrum  from 
particular  physical  models  of  IMG  sources  as  well  as  to  develop  procedures  for  diagnosing  sources  of 
IMG. 

ANALYTICAL  CONSIDERATIONS 

Intermodulation  generation  occurs  when  an  electronic  system  contains  one  or  more  elements 
which  do  not  exhibit  purely  linear  current-voltage  (i-e)  relationships.  It  is  mathematically  well-known 
and  obvious  how  non-linearities  give  rise  to  IMG. 


Consider  the  functional  relationship 


e  =  Ai  +  Bi 2 


where  e  is  the  voltage  and  /'  is  the  current  across  an  element  in  the  system. 

Then,  for  a  two-carrier  signal  with  frequencies  a>|  and  w2  and  zero  relative  phase  at  t  -  0 

i  —  /']  +  i 2  “  /]  COSttif  +  /2  cosa>2/ 


and  thus 


e  -  /4l/jCOS<oi/  +  /2cos<u2tI  +  fl[/,cos<iJit  +  /2coscu2M2 
-  A  [/,  cos«|/  +  /2  cosa>2/]  +  B[I*  cos2«if  4-  /2  cos2w2r 
+  2/|/2  COSCOjt  cos«o2/] 
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Recalling  that 

cos2*  —  1/2  (cos  2x  +  1) 

cos  *  cos  y  =  l/2[cos(x  +  y)  +  cos  (*  -  y)l  (3) 

we  find 

e  —  A  Ix  cosaj,/  +  /2cos <o2t  +  yfi{/2  (cos2aj,f  +  1)  +  /22  (cos2o>2f  +  1) 

+  2/[/2[cos (<i>i  +  w2)f  +  cos  (a>|  —  cu2) /] }  (4) 

and  the  i 2  term  is  found  to  give  rise  to  second  order  intermodulation  signals  (intermods). 

Analogously,  if  a  third  order  term  were  also  in  Eq.  (1),  Ci 3  say,  then  additional  terms  would 
appear  in  Eq.  (2),  to  wit 

y  C{/2  [cos3to,/  +  3cos<0|  r]  +  /2  [cos3<u2/  4-  3costi»2/] 

+  3/|/2  [2cos<t>|/  +  cos(<0|  +  2o)2)t  +  cos (to j  —  2a> 2) /  +  cos(a>|  —  2to2)r] 

+  3/2/2[2cOSW2r  +  COS(2t0|  +  tt>2)/  +  008(240!  —  to2)f]} 

In  principle,  this  same  kind  of  procedure  can  be  used  to  carry  out  an  evaluation  of  the  intermodulation 
signals  for  an  arbitrary  number  of  carriers  and  for  any  order  term  or  terms  in  the  current-voltage  rela¬ 
tionship. 

What  this  implies  is  that  regardless  of  the  functional  relationship  between  voltage  and  current,  if 
voltage  can  be  expressed  as  a  power  series  in  the  current,  i.e.  a  Taylor  series,  then  one  could  apply  the 
technique  used  in  obtaining  Eq.  (4)  to  predict  the  frequencies  and  amplitudes  of  the  intermods. 


A  Taylor  series  expansion  would  treat  the  current  associated  with  each  primary  frequency  as  an 
independent  variable  so  that  the  expansion  would  appear  as 


e  -  /(/)  -  /( 0) 


9/(0)  + 

9{<0)  <,+.. 
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9/2  2 

2! 

9/, 2 

+.... 


Z' kj~ 

*-i 


/(0)  +  yj- 


A  .  9 

k-l 


/(0) + -Jj- 


V  •  9 

L  ,kJi~ 

k-l  V'k 


/( 0)  + . 
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where  /(0)  implies  evaluation  of  /(/)  where  each  of  the  /*’ s  is  zero  and  n  is  the  number  of  primary 
carrier  frequencies.  Each  term  in  the  Taylor  series  could  be  treated  in  the  manner  leading  to  Eq.  (4) 
to  determine  the  resultant  intermods.  Although  this  procedure  would  be  very  long  and  tedious,  in 
principle  it  could  be  carried  out  just  one  time  for  each  term  in  the  series  and  for  various  numbers  of 
primary  carriers  and  the  results  tabulated.  We  have  carried  out  a  number  of  calculations  along  these 
lines  which  have  proved  to  be  very  important  for  learning  about  the  relationships  between  the  non- 
linearities  and  the  intermods  they  cause  as  well  as  the  relationships  among  different  intermods.  On  the 
other  hand  the  compilation  of  the  results  proves  to  be  rather  elaborate  and  could  not  include  the  effects 
of  the  relative  phases  of  the  primary  signals  at  some  particular  time.  (Consideration  of  phase  will  be 
discussed  below.) 


More  important,  however,  is  the  fact  that  some  simple  i-e  characteristics  often  found  in  real  sys¬ 
tems  cannot  be  expressed  in  a  Taylor  series.  For  example,  the  simplest  rectifier,  whose  i-e  functional 
dependence  is 

e  —  Ai  i  >  0  (6) 

e  -  0  i  <  0 
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where  A  is  some  constant,  does  not  fulfill  the  mathematical  requirements  for  a  Taylor  series  expansion. 
It  can,  however,  be  easily  treated  using  computer  techniques  and  therein  lies  one  major  advantage  of 
the  computer  analysis  as  compared  to  the  analytical  approach. 

A  number  of  useful  facts  for  an  analytical  diagnosis  procedure  are  already  obvious  from  these  dis¬ 
cussions.  First,  all  terms  in  the  i-e  functional  relationship  contribute  additively  to  the  intermod 
spectrum.  Second ,  the  highest  order  intermod  that  occurs  arises  from  the  highest  order  (non-zero) 

1  n  A  m 

term  that  appears  in  the  Taylor  series  expansion;  i.e.  the  term  — -  T  4  /( 0)  will  give  rise  to 

m  ■  k-\  vik 

m,h  order  intermods  and  in  general,  all  odd  or  even  order  intermods  lower  than  m  according  to  whether 
m  is  itself  odd  or  even. 

[As  an  as-de,  we  note  that  if  the  trigonometric  manipulations  analogous  to  those  leading  up  to  Eq. 
4  were  carried  out  for  wth  order  terms  in  the  Taylor  series  using  exponential  rather  than  conventional 
forms  of  the  cosine  function,  one  would  conclude  that  only  mth  order  intermods  arise  from  the  mth 
order  term  in  the  Taylor  series.  This  would  be  incorrect.  Use  of  the  exponential  functions  with  the 
taking  of  the  real  part  at  the  completion  of  the  calculation  is  not  valid  for  non-linear  problems.  The 
difficulty  arises  from  the  difference  between  the  power  of  the  real  part  of  a  complex  number  and  the 
real  part  of  the  power  of  the  same  number.  E.g.,  Re{(a  +  ib)2}  ^  {Re(o  +  i7>)}2.] 

Third,  it  is  trivial  to  demonstrate  that  any  particular  intermod,  say  cos(2o>,  +  <u2)t,  has  the  same 
phase  regardless  of  whether  it  arises  from  the  third,  fifth,  seventh,  etc.  order  term  in  Taylor  series. 
(We  are  assuming,  of  course,  that  the  non-linear  element(s)  giving  rise  to  the  intermods  are  found  at  a 
single  point  in  space.)  To  do  this,  imagine  the  quantity  u>xt  replaced  by  (w,  t  +  cf>),  and  similarly  for 
<o2r,  a>3t,  etc.  and  the  result  is  obvious.  On  the  other  hand,  given  a  multicarrier  signal,  it  is  likely  that 
two  or  more  different  intermods  will  have  the  same  frequency.  For  example,  if  o>j  =  10,  w2  =  11  and 
w3  =»  13,  (in  arbitrary  units)  then  2u2  -  o»|  —  o>,  +  <03  -  <o2  -  12.  In  this  situation  there  is  no  reason 
to  expect  these  two  different  intermods  to  have  the  same  or  nearly  the  same  phase.  This  provides  an 
explanation  for  an  often  seen  phenomenon  that  is  not  widely  understood;  to  wit,  a  given  intermod  is 
reduced  in  amplitude  when  an  additional  carrier  frequency  is  turned  on.  This  could  occur  for  example 
with  the  signal  whose  frequency  is  12  in  the  example  above  when  <o}  is  turned  on  if  the  a>\  +  <*>3  —  <o2 
signal  is  out  of  phase  or  nearly  out  of  phase  with  the  2 u2  -  a>\  signal. 

Fourth,  the  amplitude  or  power  (power  is  proportional  to  the  square  of  the  amplitude)  dependence 
of  an  intermod  on  the  primary  signal  input  amplitude  (or  power)  can  be  seen  from  Eq.  (5).  If  an  mlh 
intermod  arising  from  the  m,h  order  term  in  the  Taylor  series  expansion  has  a  frequency  (pjwi  +  p^u2 

+  P3“3  +  P4<*>4  +  — )  =  X  A<**/  where  the  p’s  are  positive  or  negative  integers  and  £|p,|  -  m,  then 

/  / 
the  intermod  amplitude  will  vary  as  I\PiI22l3P}l/A ...- J]  lp‘.  This  dependence  of  intermod  power 

I 

on  primary  signal  power  can  be  a  useful  tool  in  determining  the  i-e  functional  relationship.  It  must  be 
remembered,  however,  that  an  mlh  order  term  in  the  Taylor  series  expansion  will  also  generate,  in  gen¬ 
eral,  intermods  of  order  m  —  2,  m  —  4,  ....  etc.  and  all  of  these  will  also  be  of  m,h  order  in  the  fs.  For 
these  contributions,  specific  functional  dependence  of  the  amplitude  of  the  intermod  on  the  amplitudes 
of  the  various  primary  signals  cannot  be  simply  specified.  AH  that  can  be  said  without  tedious  calcula¬ 
tion  for  a  specific  situation  is  that  if  an  intermod  frequency  is  £  p,eo,,  then  the  amplitude  will  consist 

/ 

of  a  sum  of  terms  each  of  which  varies  as  J]  /,  p,+  where  £  Ip,  |  -  m  is  the  intermod  order  and  n,  is 

i  i 

an  integer  equal  to  or  greater  than  0.  Thus,  for  example,  a  three  carrier  signal  composed  of  frequencies 
oit,  <o2  and  (1)3  passing  through  a  nonlinear  device  may  well  produce  components  such  as 
cos  (2a>  j  +  <o2)t  arising  from,  say  a  5th  or  7th  order  term  in  the  current  voltage  relationship.  Although 
this  intermod  does  not  "contain”  an  w3,  it  will  have  an  amplitude  that  does  in  general  depend  on  /3.  In 
fact,  certain  experimental  observations,  such  as  the  variation  of  the  amplitude  of  a  term  like 
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cos(2<oi  +  oi2)/  with  /j,  suggests  the  presence  of  5th  or  higher  order  terms.  In  this  regard  one  should, 
however,  be  cautious.  This  could  appear  to  occur  if  a  different  third-order  term  involving  «3  happened, 
by  accident,  to  have  a  frequency  numerically  equal  to  2on  +  w2. 

Fifth ,  the  relative  magnitudes  of  the  intermods  as  a  function  of  the  number  of  frequencies 
represented  in  the  intermod,  the  number  of  carrier  signals,  and  the  order  of  the  nonlinear  term  are 
straightforward  if  tedious  to  deduce.  Suppose  the  equation  describing  the  nonlinearity  is 

e-Ai  +  Bil  +  Cis  (7) 

If  one  compares  a  three-carrier  signal  to  a  five  carrier  signal  then  the  Bi}  term  will  give  rise  to  the  same 
kind  and  magnitude  of  third  order  terms  for  both  cases.  For  five  carriers,  there  are  simply  more  third 
order  terms  arising  from  the  larger  number  of  combinations  of  frequencies  than  for  three  carriers.  On 
the  other  hand,  the  Bi 3  term  also  gives  rise  to  first  order  terms  and  these  nonlinear  generated  first  order 
terms  will  be  larger  when  there  are  five  carriers  than  when  there  are  only  three. 

Analogously,  the  C/5  term  will  generate  more  fifth  order  intermods  when  five  signals  are  con¬ 
sidered  than  when  three  are,  but  their  magnitudes  will  be  the  same  in  both  cases  for  "similai*  inter¬ 
mods,  i.e.  cos  (3co ]  —  2w2)r  is  "similar'  to  cos  (3w3  +  2o>2) t,  but  not  to  either  cos5<i>2/  or  to 
cos  Qxd  i  +  2o>2  +  <»»3 )t,  etc.  On  the  other  hand,  the  third-order  intermod  frequencies  that  arise  from 
the  Cis  term  do  have  a  magnitude  which  depends  on  the  number  of  carrier  frequencies. 

If  the  Cis  term  is  the  highest  order  term  in  the  i-e  relationship,  then  the  largest  fifth  order  terms 
will  be  the  so  called  five  carrier  fifths,  e.g.  cos(a>i  ±  w2  ±  a>3  ±  w4  ±  ws)f.  In  order  of  decreasing 
magnitude  will  come  cos(2<ui  ±  w2  ±  w3  ±  o»4)r,  cos(2o>]  ±  2oj2  ±  o>3)r,  cos  (3^,  ±  a>3)r,  cos  (3^,  ± 
a>2)r,  cos(4wi  ±  «»2)r,  cos5oj[f.  The  trend  is  clear;  the  more  frequencies  and  the  more  equally 
represented  in  the  intermod,  the  larger  that  particular  intermod  will  be. 

In  general,  higher  order  intermods  tend  to  have  smaller  amplitudes.  For  example,  the  C/5  term  of 
Eq.  (7)  will  give  rise  to  third  order  as  well  as  fifth  order  intermods,  and  in  fact  the  three  carrier  third, 
i.e.  cos  (<o |  ±  w2  ±  <o3)/  will  have  a  magnitude  somewhat  larger  than  the  five  carrier  fifth  mentioned 
above.  It  is  by  no  means  true,  however,  that  all  or  even  most  third  order  intermods  are  larger  than  all 
the  fifth  order  intermods. 

COMPUTER  DIAGNOSIS  OF  NONLINEAR  SYSTEMS 

Using  Equations  (2)  through  (4)  is  basically  a  simpler  way  of  obtaii.ing  a  Fourier  transform  than 
the  usual  analytical  procedure.  On  the  other  hand,  computer  based  Fourer  transform  analysis  is  now  a 
well  developed  technique.  This  fact  and  the  widespread  availability  of  high  speed  computers  strongly 
suggest  the  application  of  computers  to  diagnosing  nonlinear  systems  from  ti'e  system’s  intermodula¬ 
tion  generation.  The  advantages  of  using  computers  are  more  than  just  speed.  Unlike  the  mathemati¬ 
cal  requirements  on  the  current-voltage  relationships  that  are  necessary  in  the  analytical  approach  dis¬ 
cussed  above,  the  computer  approach  requires  no  particular  restrictive  mathematical  criteria  for  the 
current-voltage  dependence.  For  example,  there  is  no  special  difficulty  in  Fourier  analyzing  a  multi- 
carrier  signal  imposed  on  a  circuit  element  with  circuit  characteristics  given  by  Eq.  (6).  Thus,  given  an 
i-e  relationship,  analytical  or  otherwise,  one  can  predict  the  magnitude  and  relative  phase  of  the  inter¬ 
mods. 

There  are  three  characteristics  of  computer  based  Fourier  transforms  which  should  be  understood 
if  results  are  to' be  properly  interpreted.  First,  a  computer  based  Fourier  transform  (FT)  differs  from  an 
analytical  transform  in  that,  with  the  former,  one  works  with  a  discrete  (rather  than  continuous)  set  of 
"data"  (values  of  the  function  being  transformed)  summed  over  a  finite  (rather  than  infinite)  range  of 
the  argument  of  the  function.  As  a  consequence  of  this,  the  relative  phases  of  the  input  signals  at 
some  fixed  time  can  influence,  to  some  extent,  the  magnitudes  of  the  resultant  intermods.  This 
influence,  however,  is  quite  small. 
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Second ',  another  and  much  more  important  consideration  associated  with  signal  phase  is  that  all 
computer  results  are  numerical  and  there  is  no  way  to  distinguish  between  two  or  more  contributing 
intermods  at  a  given  frequency.  Further  and  as  mentioned  earlier,  there  is  no  reason  for  different 
intermods  of  the  same  frequency  to  have  the  same  or  any  other  particular  phase  relative  to  each  other. 
Thus  certain  characteristics  of  particular  intermods  cannot  always  be  easily  sorted  out.  This  corresponds 
to  the  actual  experimental  situation  where  various  intermods  can  overlap. 

Third,  a  possible  major  contribution  to  the  overlap  of  intermods  is  the  phenomenon  of  aliasing. 
When  N  discrete  data  points  are  separated  from  each  other  by  a  time  interval  T,  then  there  is  a  max¬ 
imum  possible  frequency,  com,  that  can  be  resolved  which  is  approximately  ~j.  for  large  N.  Frequen¬ 
cies  less  than  <am  will  be  correctly  given  by  the  computer  based  Fourier  transform.  For  those  frequen¬ 
cies  greater  than  tom  by  A<om,  the  associated  computer  generated  spectrum  line  will  appear  at  a  fre¬ 
quency  (wm  —  Aa»m)  which  is,  of  course,  within  the  range  of  allowable  frequencies.  This  is  the 
phenomenon  known  as  aliasing.  Thus,  the  probability  of  signal  overlap  in  the  range  of  allowable  fre¬ 
quencies  is  increased  simply  because  of  the  apparent  increased  density  of  intermod  signals.  This 
phenomenon  will  not  usually  present  too  great  a  problem  for  two  reasons.  First,  the  intermods  with 
frequencies  greater  than  atm  are  usually  of  higher  order  and  thus  can  be  expected  to  have  very  low 
amplitudes.  Second,  the  value  of  the  intermod  frequency  can  be  made  to  indicate  which  spectrum  lines 
have  appeared  by  aliasing.  For  example,  if  all  primary  frequencies  are  whole  numbers  but  a»maI  is,  say, 
a  whole  number  plus  0.3,  then  any  spectrum  line  centered  on  an  to  that  is  not  a  whole  number  can  be 
assumed  to  have  arisen  from  aliasing. 

In  order  to  illustrate  the  nature  of  the  intermod  response  expected  from  an  i-e  relationship  not 
amenable  to  analytical  treatment,  a  number  of  computer  experiments  have  been  carried  out.  For  this 
purpose  we  have  worked  with  linear,  quadratic  and  cubic  rectification  functions  and  various  combina¬ 
tions  thereof.  Thus,  the  i-e  relationship  is 

e  —  ai  +  bi*  +  ci3  i  >  0 

-  0  i  <  0  8 

where  any  one  or  two  of  the  coefficients  a,  b,  or  c  may  be  zero.  For  example  if  a  —  b  —  0,  c^O  then 
we  refer  to  this  as  cubic  rectification.  If  only  a  &  0,  it  is  linear  rectification  while  if  only  b  ^  0  it  is 
quadratic  rectification.  Calculations  using  up  to  three  input  frequencies  have  been  carried  out  with 

I  *“  l\  COS  (<U]f  +  <f>\)  +  1 2  COS  (o>2  t  +  <£2)  +  1 2  COS  (fiijt  +  $3) 

as  the  input  signal,  where  toi  —  10,  cu2  “  1U  and  co3  -  13  and  the  phases  <f>\,  and  4> 3  could  be  varied 
using  a  random  number  generator.  Amplitudes  of  signals  are  in  arbitrary  units  since  only  ratios  of 
amplitudes  are  significant. 

In  the  first  computer  experiment,  two  of  the  three  primary  signals  and  two  second  order  inter¬ 
mods  were  examined  for  several  combinations  of  the  coefficients  a,  b,  and  c  in  Eq.  (8).  The  parame¬ 
ters  used  and  results  obtained  are  summarized  in  Table  I.  Values  of  a,  b,  and  c  used  are  shown  above 
the  column  corresponding  to  the  results.  If  one  or  more  of  these  coefficients  are  not  indicated  above 
the  columns,  the  implication  is  that  they  are  zero.  The  results  themselves  are  averages  of  four  separate 
runs  using  random  phases  of  the  input  signals  with  each  run.  The  intent  was  to  minimize  the  effects  of 
signal  phase  discussed  above. 


Table  I 

—  Signal  amplitudes  at  various  frequencies  for 

the  kinds  of  rectification  indicated. 

Intermod  and 

C“  1 

b-1 

b~c”  1 

a—  1 

a-b— c-1 

Frequency 

wj  —  10 

34.571 

17.047 

51.618 

9.265 

60.883 

o»3  —  eu  |  “  3 

22.835 

9.153 

31.988 

3.299 

35.287 

to 2  +  a>3  ~  24 

22.697 

9.117 

31.813 

3.036 

34.848 

C03  ”  13 

34.534 

16.695 

51.229 

9.242 

60.470 
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The  most  obvious  feature  of  the  data  in  Table  1  is  the  equality  between  the  sum  of  the  amplitudes 
for  quadratic  ( b  -  1)  and  cubic  (c  =  1)  rectification  calculated  separately  and  the  result  for  the  sum  of 
quadratic  and  cubic  rectification,  column  three.  The  summability  is  also  true  for  linear  rectification. 
This  result  is  not  surprising  since  it  merely  implies  the  obvious  summability  of  spectrum  lines  from 
different  nonlinearities.  A  second  interesting  regularity  is  that  for  each  type  of  rectification,  "similar" 
intermods  have  very  nearly  the  same  amplitudes.  The  extent  to  which  the. amplitudes  are  not  identical 
results  from  more  than  one  intermod  contributing  at  a  given  frequency. 

The  variability  of  signal  amplitude  because  of  multiple  intermod  contributions  at  a  single  fre¬ 
quency,  and  the  manner  in  which  this  amplitude  can  vary  with  primary  signal  phase  are  illustrated  in 
Table  II.  Here  the  results  of  several  individual  runs  with  random  phases  are  presented  for  intermod 
frequencies  of  7  and  12.  The  reproducibility  of  the  signal  whose  frequency  is  7  is  much  better  than  the 
signal  whose  frequency  is  12  because  the  two  major  contributors  to  the  former  signal  are  third  order 
and  fifth  order  intermods.  Since  the  third  order  intermod  should  be  perfectly  reproducible  and  is  much 
larger  than  the  fifth,  whether  the  fifth  adds  or  subtracts  from  the  third  (depending  on  their  relative 
phase)  the  amplitude  of  the  net  signal  is  changed  only  to  a  limited  extent.  In  contrast,  the  major  con¬ 
tributions  to  the  signal  whose  frequency  is  12  come  from  two  third  order  intermods  whose  magnitudes 
can  be  expected  to  be  much  more  nearly  alike.  Thus,  addition  or  subtraction  (constructive  or  destruc¬ 
tive  interference)  of  the  two  results  in  large  variations  in  the  magnitude  of  the  net  signal. 

Table  11  —  Signal  amplitudes  for  theto  =  7  and  w  =  12 
signals  for  the  kinds  of  rectification  indicated. 


Intermod  and 


Frequency 

run 

C=1 

b=l 

b=c=l 

a=l 

a=b=c=l 

2ft)]  —  o>3  =—  7 

1 

2.007 

.786 

2.793 

.133 

2.876 

3a>2  —  2oj3  *=  7 

2 

2.030 

.7739 

2.803 

.137 

2.838 

3 

2.0412 

.7955 

2.837 

— 

2.839 

12 

1 

3.145 

1.422 

4.563 

... 

4.590 

2w2  —  <Of  ••  12 

2 

4.182 

1.813 

5.994 

.0418 

6.030 

3 

2.024 

1.023 

3.047 

— 

3.022 

The  surprising  thing  about  Table  II  is  that  the  sum  of  the  c  —  1  and  b  —  1  columns  equals  the 
b  -  c  -  1  column  for  an  intermod  frequency  of  12.  Unlike  the  results  of  Table  I,  the  <o  =  12  signal 
shown  in  Table  II  is  a  sum  of  two  different  intermods  of  comparable  magnitude.  It  is  easy  to  show  that 
under  these  conditions,  the  phase  of  the  net  signal  for  a  particular  nonlinearity  will  not  in  general  be 
equal  to  the  phase  of  a  signal  of  the  same  frequency  arising  from  a  second  nonlinearity.  (The  phases 
would  be  the  same  if  the  same  single  intermod  was  the  only  contribution  for  both  of  the  nonlinearities 
as  is  the  case  in  Table  I).  Thus  the  sum  of  columns  c  ■=  1  and  b  —  1  would  not  be  expected  to  equal 
column  c  —  b  —  1.  Suppose,  however,  that  Ac  is  the  coefficient  of  the  a>|  -  a>2  +  w3  intermod  resulting 
from  c  —  1  and  Mc  is  the  coefficient  of  the  2a>2  -  a>i  intermod  resulting  from  c  -  I  and  Ab  and  Mh  are 
the  analogous  coefficients  for  b  -  1.  Then,  if  Mc/Ac  -  MJ  Ab  it  can  be  shown  that  the  net  signal  from 
c  -  1  and  b  —  1  will  have  the  same  phase,  and  the  additivity  of  the  w  =  12  signals  in  Table  II  are  to  be 
expected. 

In  Table  III  a  large  number  of  intermod  amplitudes  are  given.  It  is  obvious  that  the  amplitudes  of 
the  intermods  fall  into  groupings  as  shown.  Furthermore  the  relative  magnitudes  correspond  quite  well 
with  what  we  would  expect  from  the  analytical  discussions  earlier  in  this  paper,  including  the  facts  that 
the  higher  the  intermod  order,  the  lower  is  its  amplitude  and  that  similar  intermods  have  similar  ampli¬ 
tudes.  Thus,  it  is  possible  to  check  the  assumption  of  the  previous  paragraph  that  Mc/A(  -  Mb/Ab.  By 
examining  the  magnitudes  of  the  intermods  +  o>2  +  or  w2  +  <o3  -  <o,  and  2u>\  4-  a»i  or  2 w2  +  W|, 
which  have  frequencies  at  which  no  other  intermod  of  comparable  order  contributes,  one  can  estimate 
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Table  III  —  Signal  amplitudes  of  the  first,  second  and  third  order 
intermods  for  the  kinds  of  rectification  indicated.  Intermods 
at  higher  order  are  indicated  (in  parenthesis)  when  they 
overlap  a  lower  order  intermod. 


itermod 

Order 

Intermod  & 
frequency 

c= 1 

b=l 

a  =  l 

1 

ui\  —  10 

34.595 

16.984 

9.209 

1 

Cl)  2  ~  11 

33.789 

16.704 

9.007 

1 

Cl)  3  =  13 

34.589 

16.728 

9.206 

2 

<t>2  —  CO)  =  1 

20.044 

9.044 

3.883 

2 

«>3  -  C02  =  2 

20.999 

8.966 

3.404 

2 

Cl)  3  —  ClJ )  —  3 

22.789 

9.236 

3.317 

2 

Cl)|  +  <i»2  =  21 

22.350 

9.028 

3.103 

2 

co|  +  C03  =  23 

21.961 

9.200 

3.399 

2 

a>2  +  &13  “  24 

22.579 

9.031 

3.107 

2 

(4) 

2co,  =  20  {, 

3o>2-<l>3  =  20  / 

11.764 

4.604 

1.745 

2 

ha  2  =  22 

9.843 

4.483 

2.219 

2 

(10) 

2o>3  =  26  I 

6t»4  —  4<d  1  =  26  / 

12.240 

4.641 

1.647 

3 

(5) 

Cl)  1  +  Cl)  j  —  Cl)  3  —  8  | 

3o)  1  —  hi)  2  =8  f 

13.235 

3.899 

.193 

3 

3 

(Of  —  0)2  +  (1)3  =  12) 

2cd2  —  ci)i  =  12  j 

9.566 

2.964 

.074 

3 

0)2  +  (1)3  —  0)i  =  14 

13.679 

4.328 

.111 

3 

o)(  +  0)2  +  <03  “  34 

13.601 

3.912 

.02 

3 

3o>i  —  30 

2.422 

.789 

.189 

3 

3 

3(i)2  =  33  1 

2o)]  +  0)3  “  33  / 

5.481 

1.416 

.... 

3 

3o)3  =  39 

2.208 

.516 

.146 

3 

(5) 

2o)j  —  0)2  —  71 

3(1)2  —  2(1)3  “  7) 

6.923 

1.767 

.089 

3 

3 

2(0 1  —  0)2  =  9) 

2o)2  —  W3  =  9/ 

4.350 

1.545 

.109 

3 

3 

2o)2  —  0)  1  *  12 

CD)  +  (1)3  —  0)2  =  12 

9.566 

2.964 

.054 

3 

2o)3  —  0)3  “  15 

6.550 

1.996 

.154 

3 

20)3  —  o))  =  16 

6.776 

1.665 

.136 

3 

2o)i  +  0)2  *■  31 

6.635 

1.688 

.198 

3 

2o>2  +  o)|  “  32 

6.822 

1.833 

.178 

3 

2o)i  +  0)3  =  33 

5.481 

1.416 

.... 

3 

2o)2  +  0)3  “  35 

6.649 

1.999 

.152 

closely  the  magnitude  of  Mc,  Ac,  Mb  and  Ab.  It  can  be  seen  that  the  proportion  MJ A c  —  MJ Ah  is 
valid  to  the  expected  accuracy.  It  appears  in  fact  to  be  valid  for  all  three  third  order  harmonics,  but  the 
relationship  would  not  be  expected  to  hold  in  general. 

In  Table  IV  a  similar  grouping  of  intermods  using  only  two  (different)  fundamental  frequencies  is 
shown.  These  data  represent  averages  of  only  two  separate  runs  but  nevertheless  result  in  better  repro¬ 
ducibility  of  intermod  amplitude  than  Table  III  because  of  fewer  possible  intermods  having  the  same 
frequency. 
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Table  IV  —  Signal  amplitudes  of  a  variety  of  intermod  orders 
for  the  kinds  of  rectification  indicated. 


Intermod 

Intermod  & 

c=  1 

b=  1 

a=  1 

order 

frequency 

1 

o»j  —  7 

20.680 

13.244 

9.191 

1 

I'¬ 

ll 

rq 

3 

20.687 

13.243 

9.191 

2 

o)  2  —  tuj  —  10 

15.950 

9.224 

4.985 

2 

<t>2  +  <»!  ”  24 

15.702 

9.081 

4.907 

2 

2<u,  =  14 

9.457 

4.557 

1.640 

2 

2o)2  =  34 

9.406 

4.533 

1.632 

3 

to  2  —  2<»>  i  —  3 

6.911 

2.656 

— 

3 

2a)2  —  o)  i  “  27 

6.875 

2.642 

... 

3 

(1)2  +  2o) \  ”31 

6.715 

2.580 

... 

3 

3<o,  -  21 

2.246 

.518 

— . 

4 

2a)2  —  2tO|  —  20 

2.280 

— 

.995 

4 

3«i  —  0)2  ”  4 

1.353 

— 

.330 

4 

<i»2  +  3o>|  =  38 

1.307 

— 

.319 

5 

2(i)2  —  3<o  i=13 

.380 

... 

5 

3(d  2  —  2m\  —  37 

— 

.377 

— 

5 

4«t>i  —  d>2  ”  11 

— 

.074 

— 

6 

3a)2  —  3oj|  =  30 

.252 

— 

.428 

A  number  of  other  regularities  are  noticed  in  Table  IV  where  results  for  higher  order  intermods 
are  shown.  A  cubic  rectification  (c  =  1)  does  not  show  fifth  order  intermods,  but  does  generate  3rd 
order  intermods  while  a  quadratic  law  does  not  show  fourth  and  sixth  order  intermods.  Linear 
rectification  appears  to  give  only  even  order  intermods  (except  for  the  fundamentals)  which  implies  that 
the  "third"  order  intermods  in  the  last  column  of  Table  III  probably  arise,  in  fact,  from  various  high 
order  even  intermods. 

The  effect  of  primary  signal  amplitude  on  the  amplitude  of  the  intermods  for  the  three  kinds  of 
rectification  discussed  above  has  also  been  investigated.  To  minimize  multiple  intermod  contributions 
to  a  single  frequency  we  have  used  just  the  two  frequencies  7  and  17.  To  simplify  the  interpretation, 
these  two  frequencies  are  always  taken  with  equal  amplitudes  and  varied  by  factors  of  two. 

The  amplitude  relationships  are  very  obvious  (see  Table  V).  For  linear  rectification,  if  the  imput 
signals  are  increased  by  a  factor  n,  then  each  intermod  is  increased  by  the  same  factor.  For  quadratic  or 
cubic  rectification,  if  the  imput  signals  are  increased  by  a  factor  n,  then  the  intermods  are  increased  by  a 
factor  n 2  or  n 3  respectively.  It  is  interesting  that  this  is  precisely  the  kind  of  behavior  that  is  obtained 
for  analytical  linear,  quadratic  and  cubic  relationships  although  the  particular  intermods  obtained  are 
not.  For  example,  a  cubic  i-e  relationship  (no  rectification)  would  not  produce  any  intermods  greater 
than  order  three  and  no  even  order  intermods  whatever. 

The  accuracy  of  the  relationships  of  the  amplitudes  discussed  in  the  previous  paragraph  is 
extremely  high  as  is  the  agreement  of  the  amplitudes  between  different  similar  intermods  (calculated, 
but  not  included  in  Table  V).  This  is  a  result  of  only  one  significant  intermod  contributing  to  each  fre¬ 
quency.  In  Table  VI  the  same  kind  of  experiment  is  carried  out  except  that  three  frequencies  are  used 
and,  incidentally,  the  input  signal  amplitudes  are  varied  somewhat  differently.  In  spite  of  the  numbers 
in  Table  VI  being  the  result  of  averaging  several  runs  with  random  phases  the  variation  from  the  ampli¬ 
tudes  expected  is  somewhat  greater  than  in  Table  V  where  a  single  run  is  shown. 
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Table  V  —  Signal  amplitudes  for  a  number  of  intermods  for  the  kinds  of 
rectification  indicated.  The  values  of  the  primary  signal  amplitudes 
used  are  above  the  columns  corresponding  to  the  results. 


Intermod 

Intermod  & 

1/4 

1/2 

1 

order 

Frequency 

2 

Linear  Rectification  (a  =  l) 

1 

oil  —  7 

2.2910 

4.5820 

9.1637 

18.327 

2 

a»2  4-  <o\  ~  24 

1.2157 

2.4295 

4.8594 

9.7181 

2 

2co,  -  14 

.4057 

.8113 

1.6213 

3.2433 

3 

2a>2  —  to)  —  27 

.... 

.... 

.... 

.... 

3 

3<u  ]  —  21 

.... 

.... 

.... 

.... 

4 

2a»2  —  2<U|  =  20 

.2510 

.4995 

.9998 

1.9994 

4 

3<0]  —  <«»2  “  4 

.0807 

.1630 

.3256 

.6514 

Quadratic  Rectification  (b  —  1) 

1 

a)|  “  7 

.8253 

3.3011 

13.2044 

52.818 

2 

&>2  +  o> i  «  24 

.5619 

2.2475 

8.9899 

35.9594 

2 

2oi,  -  14 

.2817 

1.1266 

4.5064 

18.026 

3 

2oi2  —  o»j  ”  27 

.1639 

.6555 

2.6223 

10.4893 

3 

3b),  -  21 

.0320 

.1279 

.5120 

2.0477 

4 

2 o)2  —  2o>!  “  20 

— 

.... 

.... 

.... 

4 

3o>j  —  o>2  ”  4 

.... 

.... 

.... 

.... 

Cubic  Rectification  (c=l) 

1 

&),  =■  7 

.32215 

2.5772 

20.6175 

164.940 

2 

o)2  +  o>i  =  24 

.2429 

1.9432 

15.5454 

124.363 

2 

2o>,  =  14 

.1461 

1.1689 

9.3512 

74.810 

3 

2o)2  —  o)i  =  27 

.1066 

.8530 

6.8239 

54.590 

3 

3o),  -  21 

.0347 

.2775 

2.21974 

17.758 

4 

2o)2  —  2o>i  —  20 

.0357 

.2858 

2.28647 

18.292 

4 

3o>i  —  o>2  =  4 

.0210 

.1677 

1.3415 

10.733 

Table  VI  —  Amplitudes  as  in  Table  V  for  three  inputs. 


Intermod 

order 

Intermod  ... 

Type 

2/3 

1 

Linear  Rectification 

1 

&>!  etc  4.59 

6.1 

9.26 

2 

wi  +  <o  2  1.6 

2.25 

3.07 

3 

<ui  +  o>2  4"  (o 3  .01 

.05 

.06 

3 

2b)  3  —  b)|  .1 

.1 

.2 

Quadratic  Rectification 

1 

b),  4.20 

7.45 

16.9 

2 

Co  i  +  ci)  2  2.25 

4.05 

9.0 

3 

b)j  +  u>2  +  W3  .42 

.75 

1.65 

Cubic  Rectification 

1 

b)|  4.32 

10.2 

34.6 

2 

b)|  4-  b>2  2.75 

6.65 

22.2 

3 

CDi  +  b)j  +  W3  1.70 

4.05 

13.6 

3 

n 

00 

3 

1 

A 

1.95 

6.85 
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An  additional  experiment  was  carried  out  to  determine  the  influence  of  the  coefficients,  a,  b  and  c 
on  the  intermod  amplitude.  One  coefficient  was  individually  varied  by  a  factor  of  two,  while  the  others 
were  held  equal  to  zero.  It  was  expected  that  the  resulting  intermods  would  vary  directly  with  the  mag¬ 
nitude  of  the  coefficient  and  this  is  exactly  what  was  found. 

CONCLUSION 

We  have  shown  here  two  approaches,  the  analytical  and  numerical,  for  obtaining  the  intermod 
spectrum  given  a  specific  nonlinearity  and  given  the  number  of  primary  frequencies  involved. 
Although  one  can  do  much  with  the  analytical  approach,  we  have  shown  that  the  computer  method  is 
the  only  way  to  handle  the  problem  for  certain  rather  common  nonlinear  systems  Nevertheless  it  is 
clear  that  the  IMG  diagnosis,  that  is  the  deducing  of  the  i-e  relationship  from  the  characteristics  of  the 
intermod  spectrum,  cannot  be  carried  out  on  a  ''prescription"  basis.  No  routine  procedure  will  quickly 
and  unambiguously  provide  the  i-e  functional  relationship  and  thus  information  leading  to  identification 
of  the  physical  origin  of  the  nonlinearity.  However,  the  results  of  the  analytical  considerations  and  the 
computer  method  set  forth  above  provide  guidance  for  procedures  that  can  be  carried  out,  on  a  case  by 
case  basis,  to  obtain  the  i-e  functional  form. 
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